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Abstract 

In this paper, we consider two-scale limits obtained with increas- 
ing homogenization periods, each period being an entire multiple of 
the previous one. We establish that, up to a measure preserving rear- 
rangement, these two-scale limits form a martingale which is bounded: 
the rearranged two-scale limits themselves converge both strongly in 

and almost everywhere when the period tends to +00. This limit, 
called the Two-Scale Shuffle limit, contains all the information present 
in all the two-scale limits in the sequence. 



1 Introduction 

Homogenization is used to study the solutions to equations when there are 
multiple scales of interest, usually a microscopic one and a macroscopic one. 
In particular, one may consider the solutions to a partial differential equa- 
tion with locally e-periodic coefficients and study their behavior as the small 
period e tends to 0. Two-scale convergence, introduced by G. Nguetseng |TU] 
and G. Allaire yj, is suited to study this particular subset of homogenization 
problems called periodic homogenization. It was later extended to the case 
of periodic surfaces by M. Neuss Radu jHl |9] and G. Allaire, A. Damlamian 
and U. Hornung [2]. It can also be used in the presence of periodic holes in 
the geometry, see [HI E] or to homogenize multilayers \V2\ \ 
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Intuitively, two-scale convergence introduces the concept of two-scale 
limit Mo which is a function of both a macroscopic variable x — also called 
slow variable — and a microscopic p-periodic variable y — also called fast 
variable — such that, in some "meaning", x H- Uo{x,x/e) is a good approxi- 
mation of Us- 

As indicated by its name, two-scale convergence captures the behavior at 
two scales: the macroscopic one and the pe-periodic one. However, two-scale 
convergence does not capture all phenomena that happens at a scale linear 
in e but only those whose length scale is pe/m where m is an integer. The 
two-scale limit of a sequence depends not only on the asymptotic scale, but 
also on the precise value of the chosen period. For example, any phenomena 
happening at the length scale of 26 will not be fully apparent in the two-scale 
limit computed with period e. The two-scale limit computed with period 2e 
will contain no less — and might actually contain more — information than 
the two-scale limit computed with period e. For example, the homogenization 
of sin(27rx/£) -|- sm{nx/e) gives a two-scale limit of uq : {x,y) sm{27iy) 
if computed with the homogenization period e, i.e., when p = 1, and Uq : 
(x, y) sin(27r|/) -|- sin(7r?/) if computed with the homogenization period 2e, 
i.e., when p = 2. Furthermore, if we choose p = 1/2, then the two-scale 
limit is none other than the null function. Worse, the scale factor p could be 
irrational. 

The choice of the scale factor p used in the homogenization process is 
therefore of utmost importance in two-scale convergence. Using a badly 
chosen scale factor p may and will often cause a huge loss of information. At 
worst, we recover no more information than the one obtained by the standard 
weak limit: if pe is the correct choice of homogenization period, the two- 
scale limit computed with period Xpe where A is an irrational number should, 
intuitively, carry no information about what happens at scale pe. 

Fortunately, there is usually a natural choice of period: the coefficients of 
the partial differential equation are often chosen locally e-periodic. The most 
natural choice is to choose p = 1, i.e., to consider the correct microscopic 
scale for is e itself. If there are two important periods to consider pe and 
p'e, the intuitive solution is to choose a period that is an entire multiple of 
both. However, this can only be done if the ratio p/p' between the two scale 
factors is a rational number. 

When the two-scale limit depends on the fast variable, we may consider an 
homogenization period of p2e instead of pie where P2/P1 is a positive integer. 
The two-scale limit computed with the homogenization period p2e contains 
more information than the two-scale limit computed with the homogenization 
period pie. It is then natural to study the behavior of the two-scale limit as 
the scale factor tends to +00. G. Allaire and C. Conca studied in [3J a similar 
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problem and established, for an elliptic problem, the behavior of the spectra 
of the equation satisfied by the two-scale limit as the scale factor p goes to 
+00. G. Ben Arous and H. Owhadi [4J studied the behavior of the Brownian 
motion in a periodic potential using multiscale homogenization when the 
ratio between two successive scales is bounded from above and below. 

In this paper, we consider various two-scale limits, each computed with 
a different homogenization period. In particular, we consider a sequence of 
periods (p„)„eN such that for all integers n, Pn+i/Pn is a positive integer and 
we study the two-scale limit of {u£)s>o computed with the homogenization 
period Pn£- This two-scale limit, denoted Uo,p„, is p„-periodic in each com- 
ponent of its fast variable. Since Pn+i is always an entire multiple of one 
can always recover the two-scale limit ito,p„ from the two scale limit Mo,p„+i- 
If Pn+i = rrinPn ^ud in dimension d > 1: 



The sequence of two-scale limits (Mo,p„)nGN yields increasing information 
on the asymptotic behavior of (Me)e>o- A natural question is whether the 
two-scale limits Mo,p„ themselves converge whenever n tends to +00. I.E., 
does there exist a function that carry the information of all the p^-two-scale 
limits? The goal of our paper is to answer this question. The answer is 
positive. We show in this paper that the sequence of two-scale limits is, 
after a measure preserving rearrangement, a bounded martingale in and 
therefore converges both strongly in and almost everywhere to a function 
we call the Two-Scale Shuffle limit. 

In ^ we remind the reader of previously known results: two-scale con- 
vergence and the convergence properties of bounded martingales. In §21 
we show how the different two-scale limits are related to each other through 
martingale-like equalities and explain how to transform these two-scale limits 
to get a genuine martingale. This leads to our stating of our main theorem: 
Theorem 13.81 in which we show that in a certain meaning the two-scale limits 
themselves converge to the Two-Scale Shuffle limit. In addition, we also state 
in Corollary 13.91 that all the information present in all the two-scale limits 
is contained in the Two-Scale Shuffle limit. In §11 we use this result on the 
heat equation in multilayers with transmission conditions between adjacent 
layers and establish, for this particular example, the equation satisfied by the 
Two-Scale Shuffle limit in Theorem 14.11 
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2 Notations, prerequisites and known results 



Throughout this paper, if x is in M, we denote by \_x\ the integer part of 
X. We also denote by |?T,i,n2] the set [ni,?7,2] H N. To make the present 
paper as self-contained as possible, we recall in this section known results on 
the two main mathematical tools we use to prove our main theorem: two- 
scale convergence in §2.1^ and classical results on the convergence of bounded 
martingales in §2.21 



2.1 The classical notion of two-scale convergence 

First, as in ^j, we introduce some notations. In this paper, p always refer to 
a scale factor. It remains constant while taking the two-scale limit. However, 
the goal of this paper is to observe the behavior of the two-scale limits as p 
tends to +oo. 

By Q, we denote a bounded open domain of M"' where d > 1. By Yp, 
we denote the cube [0,p]'^. By L'^(Yp), we denote the space of measurable 
functions defined over R'', that are p-periodic in each variable and that are 
square integrable over Yp. By C#{Yp), we denote the set of continuous func- 
tions defined on that are p-periodic in each variable. 

We reproduce the now classical definition of two-scale convergence found 
in jll [10]. For convenience, we added the scale factor p. 

Definition 2.1 (Two-scale convergence) . Let p be a positive real. A sequence 
{ue)e>o belonging to L^(i7) is said to p-two-scale converge if there exists u^^p 
in L^(f2 X Yp) such that: 

Urn I Ue{x)ij (x,-) dx = ^ I I uo,p{x,y)ij{x,y) dydx, (2.1) 
^ P JnJYp 

for all tfj in V {Vl; C#{Yp)) . 

It is a common abuse of notation to also designate by Uo,p the unique 
extension of Mo,p to x M'^ that is p-periodic in the last d variables. 

G. Allaire, see [T], and G.Nguetseng, see |1U) . proved that any sequence 
of functions bounded in has a subsequence that two-scale converges. Let's 
reproduce this precise compactness result. 

Theorem 2.2. Let {u£)s>o be a sequence of functions bounded in L^(fi). 
Then, there exist Uo,p in L^{Qx]0,p[^) and a subsequence Sk converging to 
such that 

lim I Ue,{x)ij (x, — ) dx = ^ I [ uo,p{x,y)ij {x,y) dy dx, (2.2) 
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for alltfj m V {Vl; C#{Yp)) . 

Proof. See G. Allaire [li Theorem 1.2] and G. Nguetseng [10, Theorem 2]. 
The presence of the scale factor p has no impact on the proof. □ 

We also have the classical proposition 

Proposition 2.3. Let p-two-scale converges to u^^p. Then, 

1 



Mo,p||L2(Qxyp) < liminf ||Me||L2(Q). 



Proof. See G. Allaire jl] Proposition 1.6]. The presence of the scale factor p 
has no impact on the proof. □ 



The next proposition is easy to derive from Theorem 12.21 

Proposition 2.4. Let {pn)nen be an increasing sequence of positive real num- 
bers. Let (Me)e>o be a sequence of functions bounded in L^(fi). Then, there 
exist a subsequence {ek)km converging to 0, and a sequence of functions u^^p^ 
in L'^{Qx]0,Pn[^) such that, for any non-negative integer n, the sequence 
(M£j,)fcgN Pn-two-scale converges to Mo,p„- I-E., such that for all integers n: 

lim / Ue,{x)i! [x, — ] dx = ^ [ [ uo,p„ix,y)i: {x,y) dy dx, 
for all^ m L\n;C#{YpJ) . 



Proof. Apply Theorem 12.21 multiple times and proceed via diagonal extrac- 
tion. □ 

Our goal in this paper is to study the limit of Mo.p„ as p„ tends to +oo. 



2.2 Convergence of bounded martingales 

In this section, we recall the notions of probability theory needed to prove 
our main theorem. In particular, we are interested in using the convergence 
properties of bounded martingales. For more details, the reader may con- 
sult 13 . We assume the reader to be familiar with the notions of a-field and 
cr-additivity in measure theory. 

We use the following common notations: 

• If C is a subset of V{X), we denote by a(C) the smallest cr-field in X 
that contains C. 
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• li D is a topological space, we denote by B{D) the set of all Borel sets 
in D, i.e., the smallest a-field containing all the open subsets of D. 

Definition 2.5 (Measurable space). A pair {X, J-") is said to be a measurable 
space if J-" is a a- field in X. 

Definition 2.6 (Measure space). A triplet (X, J-", /i) is said to be a measure 
space if (X, J^) is a measurable space and if is a positive cr-additive measure 



A measure space {X,J^,fi) is said to be finite if fi{X) < +00. A mea- 
sure space {X, J-", fi) is said to be a-finite if X is the countable union of 
J-'-measurable sets of finite measure. A measure space {X, J-", P) is said to be 
a probability space if P(X) = 1. 

We start by recalling the definition of conditional expectation, see jTJ 
ch. 6, Theorem 6.1] for more details. Usually, the conditional expectation 
is defined for probability spaces. The definition extends without problem to 
finite measure spaces and even, to some extent, to a-finite measure spaces. 

Definition 2.7 (Conditional expectation). Let (X, J-", /i) be a measure space 
with fi being positive and cr-additive. Let ^ be a a-field such that Q d T and 
(X, ^, /i) is also (T- finite. Let / : X — > M be J-'-measurable and in ^[^^^(X, /i). 
The conditional expectation of / with respect to the a-field Q is denoted by 
E(/|J-'), and is defined as the unique, up to a modification on a set of null 
measure, ^-measurable function g such that 



for all B in Q. 

The existence of the conditional expectation is given by Radon-Nikodym 
theorem. The measure /i need not be a probability measure. However, to 
apply Radon-Nykodim theorem, (X, ^,/i) needs to be a-finite, hence the 
restriction in the definition. A statement and a proof of the Radon-Nikodym 
theorem can be found in [11, Theorem 6.10]. 

It is not enough that (X, J-", /i) be cr-finite in Definition 12.71 

Remark 2.8. When ^ C J-", it does not follow from (X, J-", yu) being a-finite 
that (X, fi) is also u-finite. A counter-example is easily obtained by setting 
Q := {0,X} whenever /i(X) = +00. 

In our main theorem, we restrict ourselves to the case of finite measures. 
However, Remark l2.8l will explain why the martingale approach doesn't quite 



on (X,J^). 
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work for the most natural attempt to define a convergence for two-scale limits, 
see mi 

In order to define martingales, we remind the reader of the definition of 
filtration. We limit ourselves to filtrations indexed by the set N. See [7^ ch. 7, 
p. 120] for more details. 

Definition 2.9 (Filtrations). Let (X, J-") be a measurable space. A sequence 
{J^n)nGN of cT-fields, J^n C J"" is a filtration if, for all non-negative integers n, 
Tn is a subset of J-'n+i- 

We now recall the definition of martingales. 

Definition 2.10 (Martingales). Let {X,J^,fi) be a cr-finite measure space. 
Let (J-'n)neN be a filtration on {X,J-',fi) such that (X, J-o,/i) is a-finite. 

A sequence (/n)neN is said to be a (J-'„)„gN-niartingale, if for all non- 
negative integers n and j, 



We now reproduce the convergence results of bounded martingales: 

Theorem 2.11 (Convergence of bounded martingales). Let {X,J^,fi) be a 
measure space with finite measure. Let {J^njnm be a filtration on the measur- 
able space (X, J-"). Let q be in ]l,-|-oo[. Let (/n)neN be a {J^n)nm-''^o,''^'tingale 
such that the sequence {fn)n&i is bounded in L'^{X). Then, the sequence 
(/n)ngN converges both almost everywhere and strongly in L'^{X,F). 

Proof. See ^ Corollary 7.22] for the strong convergence. The almost 
everywhere convergence is stated in [3 Theorem 7.18] and holds even for 
q = 1. While these two results are stated for probability measures, the 
finite measures case is easily deduced from the probability measure case by 
considering the probability measure /i(-)/yu(X). □ 

The above theorem extends, at least partially, to cr-finite measures: 

Remark 2.12. In Theorem 12. IH if the probability space (X, J-", P) is replaced 
with cr-finite measure space (X, J-", /i) such that (X, J-'o, /u) is also cr-finite, then 
the bounded martingales converge almost everywhere and at least in L^^^. It 
is unknown to the author if the strong L*^ convergence can be generalized to 
the cr-finite case. 




/.£"., if fn is J-n-measurable and if for all F in J-'„: 




(2.3) 
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3 Two-scale limits and bounded martingales 



In this section, we always assume both of the following assumptions are 
satisfied: 

Assumption 3.1 (Integer scale ratios). We are given a real sequence (p„)„gN, 
such that for all n in N, p„ > and Pn+i is an entire multiple of p„. Moreover, 
we set for n > 1, rrin := Pn/Pn-i £ N, and for n > 0, M„ := Pn/Po G N. 

Assumption 3.2. We are given a sequence of functions (ue)e>o bounded in 
L^(r2) and a decreasing sequence of positive {sk)ken such that the sequence 
{u£f,)km ]9„-two-scale converges for all integers n to a function Mo,p„ that 
belongs to L2(r]x]0,p„['^). 

This last assumption is justified by Proposition 12.41 

Our goal is to study the convergence of the two-scale limits Uo,p„ when n 
goes to infinity. In this section, we proceed as follows: we begin by estab- 
lishing a useful equality that looks like a martingale equality in §3.1^ then 
we propose a rearrangement of the two-scale limits in §3.2[ and finally pro- 
pose another rearrangement of the two-scale limits in §3.3[ the shuffle, which 
transform the sequence of two-scale limits into a bounded martingale. 

3.1 An almost martingale equality 

We start with a simple but essential proposition. 

Proposition 3.3. Suppose both assumptions \3.1\ and \3.2\ are satisfied. Then, 
for all j in N, all n in N, almost all x in Q and almost all y inYp: 

Uo,pn{x,y) = i.^^] t^o,Pn+j{x,y + ctpn). (3.1) 

Proof. Let belong to C°°{Q x W^) be p„-periodic in the last d variables. 
Since Pn+j/pn is an integer, is also p„+j-periodic in the last d variables. 
We take the limit of J^Ui.{x)(j){x,x/e) dx, as e tends to 0, in the sense of 
two-scale convergence for both scale factors p„+j and p„: 

^11 Uo^Pnix,y)(f){x,y)dydx = 

Pn Jn JYp^ 

= I I Uo,Pn+,{x,y)4){x,y)dydx,= 

Pn+j Jn JYj,^^^ 

( X] uo^p^^^{x,y^OLp)^(p{x,y)dydx. □ 



Pt+j Jn JVf 
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The most natural approach is to consider the Mo,p„ as functions defined 
over X M'^ and to study their convergence in some meaning in 1] x W^. 
Such a convergence result would be ideal as the intuitive meaning of the 
limit would be easy to grasp. Equality (13. ip is similar to the martingale 
defining equality (12.31) . Would it be possible to use the classical convergence 
properties of martingales, see Theorem 12. 11^ to prove the existence of a limit 
to the Mo,p„? Unfortunately, the martingale approach doesn't work in this 
setting but the attempt is, nevertheless, instructive. First, we try to construct 
a filtration (J-'n)nGN for the Uo,p„- For all positive integer n, the Uo,p„ are p„- 
periodic with respect to the last d variables. Let be the set of all Borel 
subsets of X M"^ that are invariant by translation of ±p„ along any of the 
last d directions of x M.'^. Clearly, Uo,p„ is J-'„-measurable. However, the 
measure space {Q x Mf^, J^^, /i) where fi is the Lebesgue measure is not a- 
finite: any J-'n measurable subset of f2 x M*^ is either of null measure or of 
infinite measure. Therefore, the concept of (J-'n)nGN-niartingale is ill-defined, 
see Definition 12.101 and Remark 12.81 Should we attempt to verify whether 
the martingale defining equality (12.31) hold, we would get either +oo or on 
both sides of the equation. 

However, the martingale defining equality (12. 3p is satisfied if one replaces 
the Lebesgue integral of M.'^ by the limit of the mean over a ball as its radius 
tends to +00. I.E., we have for all F in J^n 

iR/n n\\ I I mx,y)eF}uQ,p„{x,y)dydx = 

= J'^'^ p^l / / l{(2;,y)eF}Mo,p„ (a;,i/)di/da;, 

R-^+co\B{Q,R)\ J^JB^^a,R) 

where -8(0, R) is the open ball of centered on and of radius R and 
where \A\is the Lebesgue measure of set A. Unfortunately, we were unable 
to derive a direct convergence result using this pseudo- martingale equality. 
To proceed further, we need to transform the two-scale limits ito,p„ in order 
to get genuine martingales. 

3.2 Rearrangement of the two-scale limits with integers 

In the previous section, we established a "martingale-like" equality for the 
two-scale limits Mo,p„- To get genuine martingales in the sense of Defini- 
tion [5]Tni we need to rearrange the Mo,p„- While we are unable to prove a 
convergence for the rearrangement of the two-scale limits presented in this 
section, the ideas behind this rearrangement provide insight on the next sec- 
tion where we introduce another rearrangement and prove its convergence. 



9 



In this section, we rearrange the tio,p„ by introducing a new variable a 
that belongs to if' . The rearrangement, denoted by fp^, depends on the slow 
variable a? G on a fast variable y G [0,po['^5 and on the new variable ol. To 
rearrange the Mo,p„ into the Wp^, we subdivide f2 x [0,pn['^ into = [pn/voY 
sets Vt X rii=i[^i''^i +Po[- Each of these sets is the product of Vl with an 
hypercube indexed by ck = (ai, . . . , a^) and we define a, ■) as taking in 

Vt X [0,po["' the same values Uo,p„ does in x ni=i['^«5 +Po[- The variable 
y represents the position of the fast variable inside each hypercube. I.E., we 
set: 



id 

{x,a,y) ^ uo,p„{x,y + poa). 
We have the following proposition 

Proposition 3.4. For all n in N, for almost all x in Q and y in Yp^, the 
cx-indexed sequence {Vp^{x,a,y))^^^d is M^-periodic in each direction of ex. 
Moreover: 

VpAx,a,y) = Vp^^^{x,a + MnP,y), (3.2) 



/3G[0,Af„+,/A/„-l]d 



for all Oi in^L^ . 



Proof. This is a direct consequence of Proposition 13.31 □ 

This in turn should encourage us to look at the following problem. 

Problem 3.5. Let's call "imbricated {Mn)n-pGi^iodic d- dimensional sequences' 
sequences that satisfy the following properties (tn,a)neN,c»ezd such that 

• for all n in N, the a-indexed sequence {tn,a)a&<i Mn-periodic in each 
direction of ex., i.e. such that for all n in N, for all ex in Z'^, and for all 
(3 in U^: 

tn,ct 'I"n,a+Mnl3) 



for all n inN, and for all ex in Z' 



tn,a. — [ ; I tn+j,a+Mnf3- 

^^^^ /3eIO,M„+j/M„-l]d 
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Study the convergence of (in,a)nGN,c»eZ'* ^ tends to +00. Under which 
condition does there exist a sequence too, a such that for aU non-negative 
integers n 

- y 

/3e[0,Af-l]'* 

or such that ^ 
or both? 



By Proposition 13. 4[ for almost all x in and in l^^, the a-indexed 
sequences (fp,^(a;, a, y))^,^^^ are imbricated (M„)„-periodic (i-dimensional 
sequences. Solving Problem 13.51 would be the first step in having a very 
elegant limit to the f p„ as a function defined on x Z"^ x Yp^. Unfortu- 
nately, we do not have an answer for Problem 13.51 While this sequence 
is morally a martingale with respect to the filtration made of the a-fields 
{a + MfiZ*^, a G |0, M„ — l]"^}, it technically is not: we have the same prob- 
lem we had in the previous section. To conclude with bounded martingales 
on the convergence, we would need a measure /i on Z"^ such that /i(Z'^) = 1, 
invariant by translation and such that /i(mZ'^) = 1/m whenever m is an 
integer different from 0. Such a measure cannot be cr-additive. If we remove 
the a additivity constraint, then /i exists: just set 



It is unknown to the author if bounded martingales converge when they are 
defined on a non a-additive measure. To avoid that problem, we introduce, 
in the next section, a different less natural rearrangement for the Mo,p„, the 
shuffle, for which we finally prove a convergence result. 



3.3 Shuffle rearrangement of two-scale limits 

In the previous section, we investigated a rearrangement where the set Q x 
[0,j9„['^ was subdivided into subsets indexed by a G Z"^. In this section, 
we finally construct a rearrangement, the shuffle, that results in a bounded 
martingale; thus establishing a convergence result for the p„-two-scale limits 
as n tends to -|-oo. To do so, we replace the variable a belonging to Z'^ 
with the variable y' that belongs to [0, If'^. Like the variable o: of the previ- 
ous section, the variable y' indicates which hypercube of edge length po we 
consider. The variable y remains unchanged and continue to represent the 
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Figure 1: One step of the measure preserving rearrangement Mi = 2 and 
M2 = 6 

location inside the hypercube indexed by y'. We set for x in Q, y in [0,po]'^ 
and y' in [0, 1]^ 

Wn{x, y, y') := Vp„{x, a{y'), y), 



where (y-(y')i = [Mny'i\ for all integers i in Using Proposition 13. 4[ we 

derive that for almost all x in Q, y in Yp^, {j, n) in N^, and cx in |0, M„ — l]'^ 

/ Wn{x,y,y')dy' = 

Jj-ld ,_ai_ aj + l , 
" lli=l[M„ ' Mn ' 



To transform the Wn into martingales, we need to shuffle the hypercubes 
as in Figure [T] where, to simplify the drawing, homogenization was only 
performed on the last component of M°', hence the presence of layers instead 
of hypercubes. In that figure, we show one step of the rearrangement. As 
seen in the drawing, each step of the rearrangement is measure preserving, 
therefore the full rearrangement is also measure preserving. We need n — 1 
such steps to fully rearrange Wn- 

To define rigorously this rearrangement, we begin by defining the function 
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Figure 2: Two steps of the measure preserving rearrangement Mi — 2, M2 
6 and M3 = 12 



that maps the rearranged layer index onto the unrearranged layer index 
RM,m -10, Mm -Ij ^ [0,Mm-l] 
i ^ M ■ {i mod m) + 



I 

m 



(3.4a) 



The application RM,m maps km + j to j'M + k when /c belongs to |0, M — 1| 
and j belongs to [0, m — 1|. We also have RM,m ° Rm,M = Rm,M ° RM,m = Id- 
Then, we set the function that maps the rearranged layer onto the unre- 
arranged one: 



RM,m{W^y'\) 



y ^ 



Mm 



+ y'- 



lMmy'\) 
Mm 



(3.4b) 



This represents only one step of the rearrangement on one component. 
For hypercubes, the permutation is the same but is done componentwise: we 
set 

W]o,i['^]o,i[", 

{y'l, ■■■,y'n)^ ih*M,miy'i), h*M,m{y'n))- 
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And obtain one step of the rearrangement on all d components. For the 
complete rearrangement on one component, see Figure |21 we set 



To get the complete rearrangement on all components we set 

i7„:[0,ir^[0,ir, 

We also have 

The function Hn shuffles the hypercubes Ylt=i[f^i/^n, {f^i + l)/^n[, hence we 
call Hn the shuffle function. 
Finally, we define 

Wn{x, y, y') := Wn{x, y, Hn{y')). (3.5) 

This measure preserving rearrangement, the shuffle, is purposefully con- 
structed so the Wn form a martingale for the following filtration of cr-fields 

J-n = B{n) X B{[0,p,]'^) X a {nti [4, ^ [,/3 e [0, Mn - If}. 

Remark 3.6. The above rearrangement of hypercubes is similar to the one 
used for computing in place the Discrete Fast Fourier Transform: the bit 
reversal. In the special case where M„ = 2", the rearrangement simply 
exchanges layers i, i.e., [i/2", (i + l)/2"[, and i', i.e., [i'/2", (i' + l)/2"[, when 
i and i' are bit reversal permutations of each other. I.E when i = Yl^Zo 

and t' = EfJo bj2^~'~' ■ 

Remark 3.7. For general M„, the rearrangement of hypercubes is also a bit 
reversal but for a mixed basis. If [M„?/'J = Yl^=i ^j-^j-i with hj in |0, mj — l\, 
then 



i-t^.^Hy' 



H:\iry:w-^+{y'-^^)\=y'- 



We now state our main result as a self contained theorem. 

Theorem 3.8 (Two-Scale Shuffle convergence). Let Q be a bounded open 
domain of with d > 1. Let {us)£>o be a bounded sequence of functions 
belonging to L'^{Q). Let (pn)nGN be an increasing sequence of positive numbers 
that satisfy Assumption \3.1i Set for alln > Mn := Pn/Po (^"^d for alln > 1 
rn^n '■= Pn/Pn-i- Lct {6k)keN be a decreasing sequence of positive real numbers 
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converging to such that the sequence {Usi^)km Pn-two-scale converges to Uq^p^ 
for all non-negative integer n. 
Set 

Wn-.nx [0,poY X [0,1]"^ ^ M 

{x, y, y') uo,p„ (a;, po [MnHn{y')\ + y) . 

where Hn is defined by Equations (13. 4p . 

Then, the sequence Wn is a bounded martingale in L^(fi x [0,po]'^ x [0, l]'^), 
in the sense of Definition \2.1l\ for the filtration 



J-„ = i3(fi)xi3([0,po]")xam 



i=l 



[0,M„-lf V. (3.6) 



And, the sequence Wn converges both strongly in L^(i7 x [0,po]'^ x [0, 1]'^) and 
almost everywhere inQx [0,pqY x [0, l]'^ to Woo, which we call the Two-Scale 
Shuffle limit. Moreover, 



Wn{x,y,y')dy' dydx = /// Woo{x,y,y')dy'dydx 





for all sets A in Tn- I-E., by Definition \2. 1\ Wn = E,{woo\J^n) ■ 

Proof. The Wn were constructed specifically so as to be a martingale for the 
filtration (13. 6p . To prove they are a martingale for the filtration {J^n)nm, we 
only need to prove that for all non-negative integer n, for almost all x in Q, 
almost all y in [0,po]°' and for all /3 in |0, M„ — l]'^, we have 



Wn+i{x,y,y')dy' = Wnix,y,y')dy'. 



lli — llMn' Mn ^ ^ ^i — l - ivin ' ivin 



I.E., we need to show that 



Wn{x, y, Hr,{y')) dy' = / Wn+i{x, y, HM„,m„+i°Hr,{y')) dy'. 



nd r Pi ft+ir JTT" \ 

d 



But Hn maps any hypercube ni=i[A/^n5 (/3i + l)/iWn[ to another hypercube 
ni=i[/5i/^n, (A' + 1)/A^n[ and Hn is measure preserving. Therefore, we only 
need to prove that for almost all x inVL, almost all y in [0,po]'^ and for all (3 
in |0,M„-lf 

/ « Wn{x,y,y')dy' = / Wn+i{x,y,HMn,m„+Ay'))^y' ■ 



lT[d tHj:±\ lr[d rj±L iiz±r 
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is satisfied. But this equality is equivalent to 



1 

^P„+i(a=,^M„,m„+i(m„+i/3 + /3'),?/) = 



which is true by Proposition 13. 4[ Therefore, the sequence Wn is a martingale 
for the filtration (J-'n)neN- 

By Proposition 12. 3[ this martingale is bounded in L^. It converges both 
strongly in L^(fi x [0,]9o]'^ x [0, l]'') and almost everywhere to a function Woo, 
see [3 Corollary 7.22]. □ 

Corollary 3.9. It is possible to recover uo^p^ from the Two-Scale Shuffle limit 
Woo- First, for all (3 m |0, M„ - if, all y' in IltitA/^n, (A + 1)/M„[, and 
almost all {x,y) in fl x [0,poY, we have 

Wn{x,y,y') = w^{x,y,y')<ly' 

because Wn = K{w^\J^n) ■ Since the shuffle function Hn is one to one from 
[0,1[°' to [0,1['^, see Remark \37!\ we have Wn{x,y,y') = Wn{x,y, H~^{y')). 
Finally, uo^p^{x,y) is equal to the constant value taken by y' i— )■ Wn{x,y — 
Po[y/Po\,y') wheny' belongs to the hypercube [[y/po\/Mn,{[y/po\+l)/Mn[. 



4 Application: heat equation in multilayers 

In this section, we consider the multilayer heat equation with three spatial 
dimensions which we homogenize along the vertical space variable, i.e., along 
the direction perpendicular to the layers. 

In [13], the author established the equations satisfied by the two-scale 
limits of the heat equation in multilayers with transmission conditions be- 
tween adjacent layers. When the magnitude of the interlayer conductivity 
between adjacent layers is weak, see [T3| §6.1], the two-scale limit depends on 
the number of layers present in the homogenization period, i.e., on the scale 
factors Pn- For given values of the slow variables {x, t), the two scale limit is 
piecewise constant in its scalar fast variable y and takes as many values as 
there are layers in a single homogenization cell. Our goal is to establish the 
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Sixth layer 



Fifth layer 



Fourth layer 



Third layer 



Second layer 
- -J-irteriayer- - 
First layer 




Figure 3: The multilayer geometry, six layers: = 6. 



equation satisfied by the limit of two-scale limits, i.e., the Two-Scale Shuffle 
limit, as defined in Theorem 13.81 

To do so, we first recall previously known results in §4.11 then derive new 
results using Two-Scale Shuffle Convergence in §4.21 

4.1 The two-scale limit of the multilayer heat equation 



We start by recalling some results we obtained in [12] ■ To avoid unnecessary 
complications, we consider here a simpler problem than the one considered 
in [T3| System (4.1)]. Let Q be Bx]0, 1[ where 5 is a convex bounded open 
subset of with smooth boundary. Let 6, < 6 < 1/2. Let / be the interval 
]6,1- 6[. For all A^, let In be [jf^^] (j + 6) /N, (j + 1 - 6) /N[. Let fi^ be the 
domain Bxl^. Let Vf^ = Bx{{j + 6)/N} and Vf'^ = Bx{{j-6)/N}. Let 

r^'+ = Uf=i' rf and r^'- = [jfJi rf ■-. Let = dBx]o, i[uro^'+ urj^-. 

Let 7 be the application on dQ^ that maps u in Hl(^]^) to its trace on 
dQ^ . Let 7' be the trace operator swapped between F^'^ and F^'": i.e., 
iu{x, {j ± 6)/N) = 'yu{x, {j =F S)/N) for all x in B and all j in |1, - 1]. 
See Figure [2] where we schematized the three dimensional domain by 
projecting it onto the two-dimensional plane. 

Let A, K and J be positive reals: A represents the heat conductivity 
inside VL^ , and J is the magnitude of the surfacic interlayer conductivity. 
For all positive integer A^, we consider the multilayer heat equation 



dt 



- A A MTV = in fi'^ X M+ (4.1a) 
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with the boundary conditions 



A 



duN 
du 







■n^Un + id' UN 

and the initial condition 



on Te X M+ 

on (r^'+ u r^'-^ 



1i7v(-,0) = U%. 



(4.1b) 



(4.1c) 



We also define the energy 

A 
2" 



K 

Vi;(aj)f da; + — / \^v\' da{x) 



B 



' N 



N 



di;. 



for all V in H^{VL^). We suppose 



SUpE^(M^) < +CX), 
TV 

the M-two-scale limit of the initial conditions tt^. We 
have the energy equality 



and denote by Uq 



E^{uj,{-,T)) + 





duN 


Jo JnN 


dt 



dx dt = E'^un{-,T)) 



for all T > 0. Because of the energy bound, for all x inQ and all j in Z/MZ, 
the function y i— )■ UQj^j{x,y) is constant in the interval ]j + 6, j + 1 — 6[. 
Moreover, that function is M-periodic in y. We denote by Uq j{x) the 
value taken by the function y Uq jvf(^) u) the interval ]j + 6, j + 1 — 6[. 

Using two-scale convergence [U |T0] and its variant on periodic surfaces [HI 
[9l [2], the properties of the two-scale limit of (u^)ivgN; solutions to the mul- 
tilayer heat system (14. ip with J ^ 0, were established in [13, Theorem 6.1]. 

Let {Mn)n>o be a sequence of positive integers such that Mq = 1 and such 
that Mn+i is always an entire multiple of For all {x,t) in Q x M"*", the 
M„-two-scale limit y Mo,Af„(iC, t, y) takes M„ values: it is constant in each 
interval ]j + 6, j + 1 — S[. For j in Z/Mn'^, we note Uo^M„,j{x, t) the value of 
Uo^Mni^jtj •) in this interval. We have Mo,M„,j+M„(a;, t) = Uo^M„j{x,t). These 
functions satisfy, for all j in TLjM^^ the weak formulation of 



dt 



2K 

A At uo,Mn,j + ^ _ 2^^o,M„,i+ 

+ YT^(2Mo,Af„,i - Mo,Af„,j+l - M0,M„,i-l) = 0, (4.2a) 
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in i7 X ]R+, and where Arp = ^ + Vt = -^V- And with boundary 
conditions 

^-^1^ = on {dBx]0, 1[) X M+, (4.2b) 

and initial condition 

«o,M„,,(-,0)=<,,„,^. ini?x]0,l[. (4.2c) 

We have recalled previously known results on the properties of the two- 
scale limits of the multilayer heat equation. In the next section, we establish 
the properties satisfied by the Two-Scale Shuffle limit of the multilayer heat 
equation. 



4.2 The two-scale shuffle limit of the multilayer heat 
equation 

We now use Theorem 13.81 to have the two-scale limits themselves converge. 
We choose M„ = 2" to avoid complications at first. We establish the follow- 
ing: 

Theorem 4.1. Letw^ be the Two-Scale Shuffle limit defined from (Mo,2")neN 
as in Theorem \3.8\ . For all {x,t) in Q x M"*", and y' in [0, 1]. The function 
w^{x,t, -,?/') is constant in the interval]S, 1 — S[. If we denote byWoo{x,t,y') 
the value ofwoo{x, t, ■, y') inside the interval ]5,1 — 5[, the Two-Scale Shuffle 
limit Woo is a weak solution to: 

(a^, t, y ) - A At Woo{x, t,y) + —Woo{x, t, y ) 

at 1 — 2d ^ A o \ 

J ^ ^ ^ ^ ' 

+Y3^(2^oo(a;, t, y') - Woo{x, t, r+(?/')) - Woo{x, t, T'{y'))) =0, 

in Q X M+x]0, 1[, and where, for all non-negative integers j : 

T+{y') = y' + 3- 2-(^+^) - 1 when 1 - 2'^ < y' < 1 - 2-'^^+^\ 
T-{y') = y' -3- 2-(^+^) + 1 when 2-(^+^) < y' < 2'^ 

with boundary conditions 
(Jill 

= on {dBx]0,l[) xm+x]0,l[, (4.3b) 

and initial condition 



Woo(-,0, ■) = w° 



oo ■ 



where is the Two-Scale Shuffle Limit of the initial conditions and where, 
as an abuse of notations, we denote by w^{x,t,y') the constant value taken 
by y w'^{x, t, y, y') in the interval ]S, 1 — 6[. 
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Proof. Let n in N and /3 belongs to |0, 2" — 1]. Let the w„ be defined as 
in §3.3^ and the Wn be defined from the 2"-two scale limits Mo,2" as in The- 
orem [321 Both the Wn and the Wn are functions defined over {Bx]0, 1[) x 
R+x]0, l[x]0,l[. For all {x,t) in (Sx]0,l[) x R+, and y' in ]0, 1[, the appli- 
cation y I—)- Wn{x, t, y, y') is constant on ]5, 1 — 5[. As an abuse of notations, 
let's also denote by Wn{x,t,y') the value of Wn{x,t,y' ,y) when y belongs to 
]6, 1 —6[. We use the same abuse of notations for the Consider a test func- 
tion belonging to C°°(nx M+). Set ^/j{x,t,y) = (p{x,t)t{y e [|r,^[}. Let 
Qt = Q X M+. Then, since the uo,2".j satisfy the weak formulation of (14. 2p . 
we have: 

Ow 

—^{x,t,y') ■ '^{x,t) dy'dxdt 
+A 1 1 I VTWnix,t,y') ■VTV{x,t) dy'dxdt 




[ Wn{x,t,y') ■ ip{x,t) dy'dxdt 
1-26 J.Iq^ %,^[ 



J 



2wn{x, t, y') ■ (p{x, t) dy' dx dt 



J 



1 — 25 hj_ £±1 

L2" ' 2" 



1 — 2(5 /[A £+lr 

12" ' 2" I 

{wnix, t, y' + 2-") + Wnix, t, y' - 2-")) ■ ip{x, t) dy' dx dt = 0, 



where, to simplify notations, we consider the function w„ to be 1-periodic in 
y'. Therefore, for all /S in |0, 2" - 1], 

" 'x, t, y') ■ ip{x, t) dy' dx dt 




+A / VTWn{x,t,y') ■VT^{x,t) dy'dxdt 

JJqtJ[4..^[ 

+ - —Ij / Wn{x,t,y') ■ ^{x,t)dy' dxdt 

+-—-T^[[ [ 2wn{x,t,y') ■ ip{x,t) dy'dxdt 

I - 2d J Jn.„ .1\A. 

l2" ' 2" L 



J 




'Qt 



1-26 J] .hjLi+i 



Wn{x, t, H*-'{H*Sy') + 2-")) ■ <^(a;, t) dy' dx dt 



J 




2" ' 2" L 

— 1/ TT*/ ^ l\ 



Wn{x,t,H:-\H:{y')-2-^))-^{x,t)dy'dxdt. = 

2^ ■'■lQTJl4r,m 



(4.4) 
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Here H* is simply the bit reversal of the first n coefficients in the binary 
expansion. Thus: 

' i 3 . 2-0+1) _ ^ <y' <1- 2-(^+^), 

r* — 1 / rr* ^ ^ ' 




K-\H:{y' + 2-")) = <( for < J < n - 1, 

if 1 - 2-" < < 1. 

(4.5) 

And 

'y' - 3 ■ 2-(J+i) + 1 if 2-(-'+i) < < 2-^ 
H:-\Hl{y') - 2-") = <; for < J < n - 1, (4.6) 

j/' + l-2-" if0<y'<2-". 

Since V9(a;,t)l{y' G ^^^[} is Jvi-measurable and Wn = '^{woo\^n), Equal- 
ity (14.41) remains valid after replacing Wn by Wao- Therefore, 



— ^(a;, t, y') ■ i^{x, t) dy' dx dt 
+A / VTWooix,t,y') ■VT^{x,t)dy'dxdt 
+ [ Wooix,t,y') ■ (p{x,t)dy' dxdt 

+ —^[ 2woo{x,t,y') ■ ip{x,t)dy' dxdt 

1 — 2o hj_ 0+l\ 

"12"' 2" L 

— ^ / {woc{x,t,T^{y')) - Woo{x,t,T'{y'))) ■ ip{x,t)dy' dxdt = 

>' Ion 1 2" L 



for all n in N and /3 in |1, 2" — 2]. Choose ?/' in ]0, 1[, for any positive integer 
n, set P = [2"?/'J and take the limit in the above equality divided by 2~" as 
n tends to +oo. □ 

If instead of setting M„ = 2", we consider a general sequence (M„)„gpj, 
the same reasonning holds. When M„ is 2", the shuffling of layers is the 
bit reversal of the first n coefficients of the binary representation of y, thus 
involutive. This is not the case for general M„ and we must use Remark 13.71 
Therefore, utmost care must be taken to compute the analogues of (14. 5 p 
and (14. 6p . We provide the limit in the general case without proof. In that 
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case, we have 



Mn 



r [y ) = y + y tt — rr — when — - — < y < — — 
for all non-negative integer j. 



K-\H:{y') + ^)={ for < J < n - 1, 



^ 1 + Yj'i=1 Ml A/|+i if < < ; 

H:-\H:{y') -j^)=\ for < J < n - 1, 

and the limit equation fl4.3ap remains valid if we set instead 



(4.7b) 



5 Conclusion 

We have proven in this paper, see Theorem 13.81 that the two-scale limits of 
a given sequence of functions, computed for periods that are entire multiple 
of the previous ones, form a bounded martingale and thus converge both 
strongly in and almost everywhere. From the limit, called the Two-Scale 
Shuffle limit, one can recover any element in the sequence of two-scale limits: 
this limit contains all the information contained in the whole sequence of 
two-scale limits, see Corollary 13.91 For a good choice of increasing periods, 
this limits captures everything that happens at any length scale that is an 
entire multiple of e. 

Unfortunately, this limit does not capture all phenomena with a period 
linear in e: it cannot capture phenomena with an irrational scale factor. 
The construction of the martingale depends on the assumption that Pn+i is 
always an entire multiple of p„. If there are two interesting scales whose ratio 
is irrational then no choice of periodic scale carry the information for both 
scales. 

We applied the notion of the Two-Scale Shuffle limit to the heat equation 
on multilayers with transmission conditions between adjacent layers. We 
then considered the solutions to these equations and established the equation 
satisfied by their Two-Scale Shuffle limit. 
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To establish the convergence of the two-scale limit, we used the shuffle 
of hypercubes described in §3.3[ Unfortunately, because of this shuffle, it is 
not easy to reach an intuitive understanding of the Two-Scale Shuffle limit. 
Results on the existence of the limit in the setting of §3.21 would not have 
that drawback. Solving Problem 13.51 would be a first step to obtain a limit 
in this setting. 
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